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Dependence of the superfluidity criterion on the 
resonance between one-particle (Bogoliubov) and 

two-particle series 



cn ■ V. P. Maslov 



'"q . Abstract 
D 

We study how the superfluidity depends on the resonance between one- and 
two-particle series. The frequency of the spectrum of two-particle solutions in an 
d . interval is calculated. 

CO 



We consider a system of A identical bosons of mass m that lie on a three-dimensional 
torus T with the side lengths L 1; L 2 , and L 2 . We assume that the bosons mutually 
\^ \ interact and that the interaction potential has the form 

§: v(nW{x-v)), (i) 

l o j ; 

where V(£) is an even smooth function rapidly decreasing at infinity and x, y are the boson 
coordinates on T. The interaction potential (JJJ depends on A such that the interaction 
radius decreases as the particle number A increases, and the average number of particles 
r**- with which a single particle interacts remains constant in this case. 

In the ultrasecond quantization over pairs, an ultrasecond-quantized operator H, 
whose explicit form was previously presented, corresponds to the boson system in ques- 
q \ tion. As stated, this ultrasecond-quantized operator satisfies identity 



H = EA, (2) 



where E is the ultrasecond-quantized unit operator and A is an operator nonuniquely 
chosen in the ultrasecond quantization space. Further, we assume that the operator A 
*P. ' has the form 



X\ A = J J dxdyB + (x, y) + ^y) + V (N l / 3 (x - y))\ B~(x, y) + 

+2 J dxdydx'dy'V (A 1/3 (a; - y)) B + (x, y)B + (x', y')B~(x, x')B~{y, y'), (3) 

where B + (x,y) and B~(x,y) are the corresponding bosonic creation and annihilation 
operators for a particle pair in the Fock ultrasecond quantization space. By identity (j2J), 
an asymptotic expression for operator (jHJ) can be found for establishing an asymptotic 
representation for the spectrum of the boson system under consideration in the limit as 

N -»• oo. 

Because the function (JTJ times A in a sense converges to the Dirac function in the 
limit as A — > oo, the operator (JHJ) involves the small parameter 1/A in the second term for 
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this limiting case. This means that to find the asymptotic expressions for eigenf unctions 
and eigenvalues of the operator A, we can use the semiclassical methods developed in 
[E] . The asymptotic expressions for eigenfunctions and eigenvalues are determined by the 
symbol of the operator Q , and this symbol is called the true symbol for the ultrasecond- 
quantized problem. The true symbol corresponding to the operator Q is the following 
functional defined for a pair of functions y) and $(x,y): 

H [$+(■), $(■)] = J J dxdy$+(x, y) (~{^> + \)j v) + 

+2 / dxdydx'dy' (NV (N 1/3 (x - y))) y)$ + (x', y')$(x, x')$(y, y'). (4) 



The conservation of the number of particles in the system for the functions $ + (x, y) and 
y) implies the condition 

dxdy® + (x,y)$(x,y) = i. (5) 

According to the asymptotic methods in j^j, in the limit as N — > oo, an asymptotic 
series of eigenfunctions and eigenvalues of the operator corresponds to each solution 
of the equations 

that also satisfy condition It follows from the explicit form of the true symbol (jlj) 
that system (jSJ can be written as 

h 2 

Q$(x, y) = - — (A x + A y )$(x, y) + 

+ JJ dx'dy' (NV (N l/3 (x - y)) + NV (N l/3 (x' - y'))) j/)$(x, x')$(y, y'), 

h 2 

n$ + (x,y) = - — (A,. + A,)$ + (x,y) + 

+ // da/dy 7 (JW (iV 1/3 (a; - x')) + NV (N 1/3 (y - y'))) $(x', y')$ + (x, x')$ + (y, y'\7) 



In the limit as iV — > oo, system (JJJ) supplemented with the condition (jSJ) has the family 
of solutions 

®t ( x , v) = t~~T2 cos ( k ( x - v)) > 

®k(x,y) = — ^ Y^^exp(zZ(x-y)), (8) 
LlL 2 j 

where and / are three-dimensional vectors of the form 



L\ L2 L12 1 

the numbers rij, n 2 , and n 3 are integers, the functions / in formula (J2J) have the form 



with the plus sign for I 2 > k 2 and the minus sign for I 2 < k 2 , and Vq here and hereafter 
denotes the expression 
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V = J dxV(x), (10) 

in which the integral is taken over the space R 3 . The vector k in (JHJ) plays the role of a 
parameter indexing different solutions of system (J2J) with conditions (jSj). The solutions (jHJ) 
are standing waves associated with the series without fluidity. 

The leading asymptotic term for the eigenvalues in the series corresponding to solu- 
tion (|HJ) is equal to the value of the product of the symbol (JIJ) on the functions (JHJ) times N, 

\2m 2 ) K J 

The asymptotic expressions for the eigenvalues and eigenf unctions, in particular, the 
asymptotic terms following (jllj) are determined by the solutions of system (J7J) and also 
by the solutions of the variational system of equations corresponding to (J7J). The detailed 
study of solutions of the variational system is here omitted. In what follows, we present 
the results of this study obtained, in particular, by using the Maple program. 

If k — 0, then the asymptotic series related to this solution (jHJ) of system (J7J), (jSJ) is the 
Bogoliubov series corresponding to the ground state without fluidity. The quasiparticle 
spectrum for this series is expressed by the well-known formula 



Xl = \l^ + V °) (12) 

If k 0, then the variational system equations corresponding to (J7J) gives the quasi- 
particle spectrum 

h 2 



= ± 



12 ]2 i / / AmV 

— \p 4 + 2 + ! - m2 - m * + # + z ? - 2jfc2 ) v (/2 ~ 11)2 + { 



h 2 



h 2 PI 2 1 / /4-mVn\ 

±— * k* + - + 1 - m 2 - m\ - i(p + q - 2k 2 ) ] j(i 2 - 1 2 ) 2 + l^-^j (,13) 

where l\ = I + 2k and I ^ k. In the formula for Ai,*^ the plus sign is chosen for I 2 > k 2 
and the minus sign for I 2 < k 2 ; in the formula for \2,k,i, the plus sign is chosen for l\ > k 2 
and the minus sign for l\ < k 2 . 

The explicit form (|13|) implies that this quasiparticle spectrum contains some negative 
values, which means that the series corresponding to the solution (JHJ with k ^ is non- 
metastable. In this connection, we suggest the following explanation for the dependence 
of the critical velocity on capillary width. In what follows, we assume that L\ 3> L 2 . We 
consider the Bogoliubov series corresponding to the flow with the velocity hk Q /m along 
the capillary, where k = 27r(n 1 /L 1 , 0, 0). For the boson system in question, the leading 
asymptotic term for the eigenvalues in this series is 

Nm + ^). d4) 



We now assume that there is a relation between L\ and L 2 such that there exists a vector 
k = 27r(0, 112/ L2, n^jL-i) for which the corresponding value of the term in (jlljl is exactly 
equal to Q14j) . This means that there can be resonance between the current states of the 
Bogoliubov series and the states of the nonmetastable series corresponding to (|S|). If the 
value of Li is very large, then resonance is also possible for the case where the value of 
the term in is close to expression ()14|) but need not coincide with it. The existence 
of such a resonance indicates the possibility of a transition from the current state to a 
nonmetastable state from which the system falls to the lowest energy level, which indicates 
the loss of superfluidity. 

The following versions of series of two-particle solutions are determined by the relation: 

Q(ki, k 2 , mi, m 2 , m 3 , h, h, h) = —5-. 2 , 3 , 2 , f 2 — s — n > ( 15 ) 

m\ + m 2 + k{ + k 2 — /f — / 2 

where fci, fc 2 , /1, / 2 , / 3 , m 1; m 2 , m 3 are integers. We impose the following conditions on the 
numbers ki, k 2 , h, / 2 , I3, m 1; m 2 , m 3 : 

Q > 0, fc? + k\ > 0, (16) 

where Q is an irreducible fraction, i.e., its numerator and denominator do not have com- 
mon divisors. 

Professor A. A. Karatsuba kindly calculated for me the frequency of the spectrum of 
two-particle solutions in the interval [a, a + h]. The asymptotic expression in the limit as 
L — » 00 has the following form. 

We are given the functions ()15j) and the function 




V = \% + ll+%, (17) 



where ki, k 2 , mi, m 2 , m 3 , li, l 2 , 13 are nonnegative integers that do not exceed L; moreover, 
Q > and k\ + k\> 0. 

We are also given the numbers a > 1, < h < 1. 

Let K be the number of sets ki,k 2 ,mi, m 2 m 3 , l\,l 2 , h satisfying the inequality 

a<V<a + h. (18) 
Then K is determined by the approximate formula 

K - i\ + 5) m L2{{a + ft > 6 " a,,) ^ ls i2( < a + ft > 6 " ' /y> ' < 19) 



If the distance between the spectrum points is of the order of jr, i.e., if it is determined 
with the desired accuracy, then the Bogoliubov series is in resonance. 

This readily shows that a ~ 4, and hence resonance occurs approximately at the points 
8nh/mL 2 . This resonance is braking the fluid, and the superfluidity disappears. 

We here present A. A. Karatsuba's proof of formula (|19j). 

We introduce the new notation 

r = ml + m 2 2 + k\ + k\, I = l{ + l\, £ = ^ > 0. 

'3 



In this notation, we have 



1 - f 2 

Q = ll -r > 0. 

A r - I 
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Two cases, (I) and (II), are possible: 

(I) < £ < 1, I < r; (II) £ > 1, / > r. 
The function V can be rewritten as 



and hence becomes 

a> <V 2 = r -f^<(a + h) 2 . (20) 

Case (I). From (j20j) we obtain 

a 2 (l - £ 2 ) + /£ 2 < r < (a + h) 2 (l - £ 2 ) + /£ 2 , 

l<r, 0<£<1. (21) 

If in (|21j) we have the relation 

a 2 (l -£ 2 ) + Z£ 2 < Z, i.e., a 2 <l, 

then system (|21|) is equivalent to 

l<r < (a + /i) 2 (l-£ 2 ) + Z£ 2 , 
a 2 <l < (a + h) 2 , 

< £ < 1. (22) 

But if in (12 1|) we have 

a 2 (l -£ 2 ) + /£ 2 > Z, i.e., /<a 2 , 

then system (}2"T]) is equivalent to 

a 2 (l-£ 2 ) + Z£ 2 <r<(a + /i)(l-£ 2 ) + /£ 2 , 
I < a 2 , 

< £ < 1. (23) 



Let K\ be the number of sets satisfying either (j22|) or 
It follows from systems (jUj) and (|23Jl that 



/ < a 2 or / < (a + /i) 2 , 

i.e., we always have l\, 1% < a + h < a + 1, since /i < 1. 
Precisely in the same way we obtain 

r < (a + /i) 2 , i.e., mi, iri2, ki, k 2 < a + h < a + 1. 

Therefore, the six variables h,h, mx,m2'mki, k 2 are bounded by the value a + 1, and they 
play an unimportant role for small a. The main parameter is L, L — > +oo, and this 
parameter is related only to the variable £, 

£ = ^ < 1, < m 3 < h < L. 



The ranges of these variables are determined automatically by the first two inequalities 
in systems (|22jl and (|2Hj) . In what follows, we present two asymptotic formulas (in the 
limit as x — > oo). 

The number of sets m\, m 2 , k\, k 2 for which r = m 2 + m\ + k\ + k\ < x is equal to 
where ^(-y/i) is the volume of a four-dimensional sphere of radius y/x, i.e., 

7T 4 / 2 7T 2 

u 4 (V5) = (v^) 4 ^4 7 = x 2 — ; (24) 

; vv ; T(4/2 + 1) 2 ' v ; 

the number of sets /1, 12 for which Z 2 + l\ < x is equal to 

— TTX. (25) 

4 

These are approximate formulas; the larger x, the more precise they are. 
From (|2I1) and ((22), 423) we obtain 

#i=£ E' ^{((«+^) 2 ( l -a+^ 2 ) 2 -(o 2 } 

0<£<la 2 <«<(a+h) 2 

+ E E / ^(( a+/i ) 2 ( 1 -^ 2 )+^) 2 -( a2 ( 1 -^)+^) 2 )}- 

0<C<1 Ka 2 

The prime on the sum over I means that the numbers I are regarded with multiplicity 
taken into account (the number of solutions of the equation is I = l\ + Z|, l\ > 0, Z 2 > 0); 
the multiplicity is equal to 7r/4 (averaged over /), i.e., 

Ki= E E r^ {((a + ^ )2(1 ~ e) + ief ~ l2} 

0<£<1 a 2 <l<(a+h) 2 

+ E T,^ 8 {{{a+h)2{1 ~ e)+ie)2 - {a2{1 - e)+ie)2)} 

0<£<ll<a 2 

The sum over I (with a good accuracy) is equal to the integral 

R i = m E / 2 (a+/l)2 {((a+^) 2 (i-a+^) 2 -^}^ 

o<£<i 

+ E r 2 {(( fl + ^ 2 ( x - « + W " - ?) + ^) 2 > dl 
o<c<i ^° 

0<£<1 v 7 



Since 



£ = < m 3 < Z 3 < L, 

'3 



the sum over £ (with a good accuracy) is equal to the integral 



Hence we obtain 
Case (II). 

£ > 1, I > r, i.e., Z 3 < m 2 , r < I. 

Let 

< £i = - = < 1, r < Z. 
S m 3 

Then we have the equality 

and the condition a < V < (a + h) implies 

a 2 < l ^$<{a + hf 
1 si 

or 

a 2 (l - ex) + Hf < I < (a + hf{l - g) + r£ 2 , 

r<l, < & < 1. (26) 

It follows from J25J) that r < (a + hf . If r > a 2 (l - £ 2 ) + rf?, i.e., if r > a 2 , then ((211) 
is replaced by the system 

r _ Z< ;<( + /l )2 (1 _ e 2 )+re 2 ) 

a 2 < r < (a + h) 2 , < & < 1. (27) 

If r'a 2 (l - £ 2 ) + r£ 2 , i.e., if r < a 2 , then (J2EJ) is replaced by 

a 2 (l - £ 2 ) + ril <l<{a + h) 2 {l- g) + r£ 2 , 

r < a 2 , < 6 < 1. (28) 

It follows from Case (I) and systems (f2*2*j) and (f2"HJ) that £ is replaced by £i, while r and Z 
interchange. 

Therefore, we use formulas (J25J) . (j2ZJ), and (J28J) to find i^, which is the number of sets 
satisfying (l2"?J) and (l2"KJ) : 

^ 2 = r E E' ((« + ^) 2 (!-^)+^ 2 -r) 

0<§i<la 2 <r<(a+/i) 2 



+ E E'^+^-^x 1 -^ 2 )- 



4 

0<£i<l r<a 2 

The prime on the sum over r means that the number of solutions r is regarded with 
multiplicity taken into account (the number of solutions of the equation is r = m\ + m\ + 
k\ + fcf). 

The number of solutions of the inequality 

m\ + m\ + k\ + k\ < x, mi, rri2, ki, ki > 0, . 
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is approximately equal to 



32* • 



We obtain 

E' - E' E 1 = E E 1 

a 2 <r<(a+/i) 2 a 2 <r<(a+/i) 2 0<t<r 0<<<(a+?i) 2 max(a 2 ,t)<r<(a+?i) 2 

= E E 1+ E E 1 

0<t<a 2 a 2 <r<(a+/i) 2 a 2 <t<(a+/i) 2 t<r<(a+h) 2 

r 2 



7T 

32 



Vn^/. H12 „2^*^/„ i 1,12 / 



■ 0<t< (a+fc) 2 a 2 <t< (a+/i) 2 



J- ( (a + /i) V - a 6 + (a + /i) 4 ((a + /i) 2 - a 2 ) - J((a + />) 6 - a 6 ) 



;((a + /i) 6 -a 6 ); 



3-16 

We derive the following formula for K 2 : 



E {(« + ^) 2 (l-ei 2 )^(0 4 -a 4 ) 
<6<i ^ 



" (1 ~ ei2) 3~16 (( ° + kf " ° 6) + (( ° + ^ " ° 2)(1 " ^32°' 

= y^((«+^) 6 -« 6 ) E (i-^ 2 )- 

0<5i<l 

Since £i = ^, < Z < m < L, we have 



o<a<i 

and 



/\ 2 \ ... L 2 
3 



The fact that K = K\ + K 2 implies the relation 

*' = (} + 5) w^ iLH{a + ,!)6 - a6) = ls i2((a + ft)6 " a% 

Remarks. 

1. The formula for K is an approximate formula. 

2. This formula is sufficiently exact if L is large and a is not too small but significantly 
less than L; for example, a = y/~L, L — > +00. 

3. If a is chosen to be not large, for example, a = 4, h = |, or a = 3, /i = |, or a = 5, 
/i = |, then we can obtain sufficiently exact formulas for L > 8, since, in this case, we 
can exactly calculate the multiplicities of / and the / themselves (as sums of two squared 
numbers). But, in this case, we must perform specific calculations. 

Example of calculations for small a. 
Let a — 4, h — 1; L > 5. 



In system (J22J), I satisfy the condition 

16 < I = l{ + l\ < 25. 
The pairs (I1J2) satisfying this condition are: 

(0,4), (1,4), (2, 4), (3, 3), (4,1), (4,1), (4,0), 



i.e., from ()22|) we obtain two systems for I = 16, I = 17, and I = 20 and one system for 
I = 18: 

16 < r < 25 - 9£ 2 , 

< f < 1, (29) 



17 <r <25-8£ 2 , 

< £ < 1, (30) 
20 < r < 25 - 5£ 2 , 

o<e<i, (31) 

18 < r < 25 - 7£ 2 , 

< £ < 1. (32) 

To calculate the number of sets (mi, 777.2, ki, k-z) for which r = m 2 + 777.3 + ^1 satisfies 
systems (|2I2J)-(JH2j) , we can either use the computer software (in this case, we must divide 
£ into intervals, where 9£ 2 lies between two neighboring integers from to 9, the same 
for 8£ 2 , 5£ 2 , and 7£ 2 , and calculate the number of sets directly) or use the formula for 
the number of integer points in the four-dimensional ball and thus find the points with 
positive coordinates. This is an approximate formula and if 

m\ + m 2 + k\ < x, rrii, ki > 0, 

then this number is equal to 

^ 2 

x 2 . (33) 



32 

Therefore, for system (|22|). we obtain the number of sets 



TV 2 



32 

0<£<l 



{((25 - 9£ 2 ) 2 - 16 2 ) • 2 + ((25 - 8£ 2 ) 2 - 17 2 ) ■ 2 

D<£<1 

+ ((25 - 5£ 2 ) 2 - 20 2 ) • 2 + ((25 - 7£ 2 ) 2 - 18 2 )}. (34) 

Since the number of such / is equal to 2, we use the factor 2. In system ()23j) . we have 

l = l\ + ll< 16, 

and the pairs (h, l 2 ) are: (M), (0^2), (M), (0^0), (1^2), (M), (^0), (2^0), (2^1), (2, 2), 
(2,3), (3,0), (3,2), (3,2). 



1= 





1 


4 


9 


5 


10 


8 


13 


multiplicity 


1 


2 


2 


2 


2 


2 


1 


2 



System ()23j) for each I has the form 

16(1 - £ 2 ) + 1? < r < 25(1 - £ 2 ) + Z£ 2 , < i < 1. 

Using (jUIH), we obtain the following formula (with multiplicity of I taken into account) for 
the number of sets for system 



7T 

32 



o<€<i 



(25(1 -a)' -(i6(i -a) 2 



+ ((25(1 -a 

+ ((25(l-£ 2 

+ ((25(1 -a 
+ ((25(l-a+5£ 
+ ((25(l-a + 10£ 



+ ((25(1 
+ ((25(1 



e 2 



a 2 - 


(16(1- 


a- 


-a 2 ) -2 


4a 2 - 


-(16(1- 


-a 


+4a 2 )-2 




-(16(1- 


-a 


+ 9£ 2 ) 2 ) ■ 2 


5£ 2 ) 2 - 


-(16(1- 


-a 


+5a 2 ) "2 


ioe 2 ) 2 


- (16(1 


-e 2 


) + 10a 2 ) • 2 




-(16(1- 


-a 


+ 8f 2 ) 2 ) • 1 


13£ 2 ) 2 


"(16(1 


-e 


) + 13£ 2 ) 2 ) • 2} 
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|(25(i -a) 2 -(16(1 -ay 

0<£<1 

+ ((25 - 24^ 2 ) 2 - (16 -15^ 2 ) 2 ) 
+ ((25 - 21^ 2 ) 2 - (16 - 12^ 2 ) 2 ) 
+ ((25 - 16£ 2 ) 2 - (16 - 7£ 2 ) 2 ) ■ 
+ ((25 - 20a 2 -(16-ll<T 
+ ((25 -15a 2 -(16 -6a 2 ; 
+ ((25-17£ 2 ) 2 -(16-8£ 2 ) 2 ; 



■-2\2\ 



• 2 

• 2 

2 

• 2 
2 

1 



+ ((25-12£ 2 ) 2 -(16-3£ 



~2\2\ 



(35) 



Since £ = < m 3 < i 3 , 1 < l 3 < L, the number of sets corresponding to Case (I) 
is either calculated directly according to (|M|) - (j33|) or is approximately replaced by an 
integral of the form 



l 

dl I dm . 



equal to 



where 



/i 




9m 



2\2 



64 P ' 



16 2 ) + 2((25 - 8m 2 ) 2 - 17 



20 2 ) + ((25 - 7m 



2\2 



u z )Y + 2((25 - 2Au z Y - (16 - 15m 
2((25 - 16m 2 ) 2 - (16 - 7m 



,2\2\ 



2((25 

+ 2((25-5m 2 ) 2 - 

+ (25(l-a) 2 - 
+ 2((25-21m 2 ) 2 

+ 2((25 - 20m 2 ) 2 

+ 1((25 - 17m 2 ) 2 - (16 - 8m 2 ) 2 ) + 2((25 - 12m 2 ) 2 - (16 - 3m 2 ) 2 )} du. 



(16(1 

- (16- 

- (16- 



12m 



2\2\ 



.2\2\ 



11m 2 ) 2 ) + 2((25- 15m 2 ) 2 - (16 



6m 



2\2\ 



in 



Thus, we have 




Case (II), i.e., the calculation of K 2 is the same; only the multiplicity of r, 16 < r < 25, 
r = m\+m 2 +k 2 +k 2 , is calculated either by an exhaustive search of the sets (mi, m 2 , k\, k 2 ) 
(also for r < 16) or by the formula for the number of integer points in the four-dimensional 
ball; the number of sets for I < r from ()27jl and ()28jl is calculated approximately by the 
formula 

7T 

r- 

Therefore, if *c(r) is the multiplicity of r, then 

K 2 = { E ^)^(25-r)(l-e i 2 ) + Er( r )( 25 - 16 )( 1 -^)} = 

0<£i<l ^ 16<r<25 r<16 ' 

= T?i( E «(r)(25-r) + 9X)«(r)). 

^16<r<25 r<16 1 

But we have already calculated that 

E K ( r ) 

r<16 

E K ( r ) 

16<r<25 

/t(r)r 

16<r<25 

Hence we have 

K 2 = -L 2 (25--(5 4 
12 V 32 v 

K = K 1 + K 2 = cL 2 . 

All these formulas are approximate. 

The number fi = j Q . . . dm and the constant c in the formula for K have not been 
calculated, but these calculations are very simple. 

The author wishes to express his deep gratitude to Professor A. A. Karatsuba and to 
G. V. Koval' for their help in preparing this paper. 
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